pss header will be provided by the publisher 



Magnetic relaxation in a model of interacting nanoparticles in 
terms of microscopic energy barriers 

Oscar Iglesias * 1 and Amflcar Labarta 1 

1 Departament de Ffsica Fonamental, Universitat de Barcelona, Diagonal 647, 08028 Barcelona, Spain 

Received May 2004, revised XXX, accepted XXX 
Published online XXX 

PACS 05.10.Ln,75.40.Mg,75.50.Tt,75.60.Lr 

Monte Carlo simulations are used to study the magnetic relaxation of a system of single domain particles 
with dipolar interactions modeled by a chain of Heisenberg classical spins. We show that the so-called 
Tln(t/ro) method can be extended to interacting systems and how, from the computed master relaxation 
curves, the effective energy barrier distributions responsible for the relaxation can be obtained. A transition 
from a quasi-logarithmic to power-law behavior of the relaxation as the interaction strength is in-creased 
is found. By direct computation of the effective energy barriers of the system, we show that this is due to 
the appearance of an increasing number of small energy barriers caused by the reduction of the anisotropy 
energy barriers as the local dipolar fields increase. 
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1 Introduction 

Dipolar interactions are ubiquitous in most magnetic nanoparticle systems due to their long-range char- 
acter. Many of the peculiar or anomalous phenomena in this kind of systems are direct consequence 
of the predominance of dipolar interactions over exchange interparticle interactions. Current preparation 
techniques and materials synthesis allow to study the role of dipolar interactions on the magnetic proper- 
ties in a controlled manner. For example, in granular metal solids the variation of volume metal fraction 
and size of the clusters allows to tune the strength of interparticle interactions Q). A similar control can 
be achieved in ferro fluids by dilution of the magnetic particles embedded in it Q. More recently, one- 
dimensional arrays of nanoelements have also been the focus of several studies |3|. In particular, there is 
current controversy on whether the slow relaxation and blocking phenomena observed in particle systems 
is to be ascribed to dipolar interactions or is indicative of a spin-glass-like state at the particle sur face |4|. 
Here we present the results of Monte Carlo (MC) simulations aiming to clarify the role played by dipolar 
interactions in the time dependent magnetic properties of nanoparticle systems, establishing a connection 
between the microscopy energy landscape of the magnetic system and the observed relaxation laws by 
means of the so-called rin(t/r ) scaling. 

2 Model and simulation method 

The model consists of a linear chain (along the x axis) of N classical Heisenberg spins S, which are 
intended to represent the a monodomain particle with magnetic moment /i; = /iSj. The particle easy-axis 
directions n; are randomly distributed and the anisotropy constants Kf are distributed according to a log- 
normal distribution of width a — 0.5 and mean value Kq = 1. Therefore, the energy of the system can be 
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H + H? ip is the effective field at site i, H is the external magnetic field and 
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is the dipolar field felt by spin Si, where g = /io/i 2 /47ra 3 characterizes the strength of the dipolar interac- 
tion. Simulations of the relaxation of the magnetization have been performed by means of the MC method 
with standard Metropolis algorithm. As for the details, let us only mention that, in order to be able to 
simulate successfully the long-time relaxation of the magnetization, we consider trial jumps between spin 
orientations corresponding to energy minima. Therefore, the energy jumps appearing in the acceptance 
probability always correspond to real energy barriers, at difference with most implementations of the MC 
method for continuous spins. In order to be able to compute the exact energy minima and barriers from 
Eq. ^ spins are restricted to point inside the x-z plane. In this way, we have been able build an algorithm 
that computes exactly the maxima and minima of the energy as well as the corresponding energy barriers. 

3 T ln(t/r ) scaling and effective energy barrier distributions 

In order to simulate the relaxation curves, we start from a configuration in which the spins are aligned in 
parallel along the z axis, a situation which mimics the application of a saturating magnetic field. Before 
starting the simulation, the spins are submitted to a previous equilibration process at T = 0, during which 
they are consecutively reoriented along the nearest energy minima during a number of MC steps until the 
system reaches an equilibrated state. Starting from this state, the relaxation towards the equilibrium state 
in zero applied field and finite T has been computed for values of g between 0.1 and 0.5 during at least 
10000 MC steps. The obtained relaxation curves have been analyzed following the Tln(t/ro) scaling 
method introduced in our previous works (5J |6| [7) for non-interacting systems. The master relaxation 
curves obtained after scaling of the curves at different T along the horizontal axis by multiplicative factors 
T are presented in Fig. ^for a range of temperature covering one order of magnitude. 

Let us first notice that the curves overlap in a wide range of times and temperatures, this overlap being 
better for the lowest T relaxations, while for the high T curves it is accomplished only above the inflection 
point. This result supports the validity of the scaling in the presence of dipolar interactions, and gives 
support to experimental observations in systems where interactions cannot be neglected 1 8 1. When dipolar 
interparticle interactions are considered, the energy barriers responsible for the thermal relaxation change 
as the relaxation proceeds but, as we have checked, this does not preclude the accomplishment of the 
Tln(t/ro) scaling. In order to understand this, we have computed, within the scope of our model, the 
evolution of the distribution of energy barriers at different stages of the relaxation process and proved 
that, in spite of the dynamic change of the local energy barriers, they remain constant during the relaxation 
process. The second point studied has been the influence of the dipolar coupling strength g on the relaxation 
law. For this purpose, we have plotted in Fig. [2] the master relaxation curves for g ranging from 0.1 to 
0.5 after smoothing and filtering of the curves in Fig. [2 Inspection of the curves shows that there is a 
change in the relaxation law with increasing g. Whereas for weak interaction (g = 0.1, 0.2) the curves 
show an inflection point around which the decay is quasi-logarithmic, in the strong interaction regime 
(g = 0.3, 0.4, 0.5), the curves have downward curvature and can be fitted to a power-law decay of the form 
m(t) oc with an exponent that decreases with increasing g, 7 = 1.02, 0.80, 0.74 for g = 0.3, 0.4, 0.5 
(dashed lines in the inset of Fig. [2] This power law behavior has also been observed experimentally in 
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Fig. 1 Left panels: relaxation curves for temperatures ranging from T — 0.02 (uppermost) to 0.2 (lowermost curves) 
in 0.02 steps for dipolar interaction strength g = 0.1, 0.4. Right panels: master curves corresponding to the relaxations 
shown in the left panels obtained after scaling by the respective multiplicative factor T. 



arrays of micromagnetic dots |9| and discontinuous multilayers 1 10 1 and also in other MC simulations 
Ising spins 1 1 1 1 and two dimensional spin systems 1 12 1. 



of 




Tln(t/x ) 



Fig. 2 Master relaxation curves for dipolar interaction strengths g = 0.1, 0.2, 0.3, 0.4, 0.5. The inset shows the same 
curves in a logarithmic scale together with the fittings to m(t) oc t" 1 (dashed lines). 



In order to better understand the origin of this change in the relaxation laws, we will proceed to extract 
some information about the microscopic energy barriers responsible for the thermal relaxation in both 
regimes. For this purpose, we will use the master curves previously obtained to apply the method proposed 
previously for non-interacting systems 0. By performing the logarithmic time derivative of the master 
curves S(t) — dM (t) / dln(t) , it is possible to infer the effective energy barrier distribution feff(E) that 
would give rise to the relaxation curve corresponding to the master curve. Fig. [5]presents the results for 
the same values of g as in Fig. [2] The main features observed are the broadening of the energy barrier 
distribution and the increasing contribution of small energy barriers as the dipolar interaction g increases. 



© 2003 WILEY- VCH Verkg GmbH & Co. KGaA, Weinhcim 



6 



6. Iglesias and A. Labarta: 



For weak interaction values (g = 0.1), f e ff(E) has a shape very similar to that for the non-interacting 
case, but its width increases with g and the mean effective barrier corresponding to the maximum shifts to- 
wards lower energies. However, when entering the strong interaction regime (marked by the appearance of 
a contribution of almost zero energy barriers), the effect of the interaction no longer consists in a distortion 
of the original distribution. In this case, the interaction creates high energy barriers resulting in a more 
uniform f e ff(E), that now becomes a decreasing function of the energy extending to higher values of 
the energy. In order to unveil the information given by the f e f f (E) in the interacting case, we have also 
computed the cumulative histograms of energy barriers that have been jumped during the relaxation process 
simulated by the MC method. Results obtained after 10000 MCS are shown for two values of g in the left 
panels of Fig. [5] and are compared to the distributions obtained from the master curves (dashed lines in 
the right panels of Fig. [3] If the histograms are computed by counting only the energy barriers that have 
not been already jumped up to a time t, we see that they coincide with the distributions deduced from the 
master relaxation curves. The differences between both results at high energy values are due to very high 
energy barriers that can only be surmounted at much longer times for the temperature considered here. 




Fig. 3 Left panel: effective energy barrier distributions obtained form the derivatives of the master curves shown in 
Fig. |2|for several values of the dipolar interaction strength g. Right panels: cumulative histograms of jumped energy 
barriers during a relaxation process at T = 0.1 after t — 10000 MC steps (symbols). Dashed lines stand for the 
derivatives of the master curves shown in the left panel. 
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